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Abstract 

We show how the early data runs of the LHC can provide valuable checks of the dif- 
ferent components of the formalism used to predict the cross sections of central exclusive 
processes. The 'soft' rapidity gap survival factor can be studied in electroweak processes, 
such as M^+gaps events, where the bare amplitude is well known. The generalized gluon 
distribution, in the appropriate kinematic region, can be probed by exclusive T produc- 
tion. The perturbative QCD effects, especially the Sudakov-like factor, can be probed 
by exclusive two- and three-jet production. We discuss the possible role of enhanced ab- 
sorptive corrections which would violate the soft-hard factorization implied in the usual 
formalism, and suggest ways that the LHC may explore their presence. 



1 Introduction 

Central exclusive production is now recognized as an important search scenario for new physics 
at the LHC, see for instance, [1] - [7] and references therein. The experimental studies of such 
processes are at the heart of the FP420 project [8, 6] which proposes to complement the CMS 
and ATLAS experiments at the LHC by installing additional forward proton detectors 420m 
away from the interaction region. In particular, these detectors will allow the measurement 
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Figure 1: A symbolic diagram for the central exclusive production of a system A. 

of the exclusive production of new heavy particles, such as Higgs bosons. As demonstrated in 
[9] - [12] such measurements will be able to provide valuable information on the Higgs sector 
of MSSM and other popular BSM scenarios. 

Indeed, central exclusive processes are very interesting both from the viewpoint of theory, 
since they contain a mixture of soft and hard QCD effects, and of experiment, as they provide a 
clean environment to measure the quantum numbers and masses of new objects which may be 
seen at the LHC. Moreover, the Jz = selection rule (J^ is the projection of the total angular 
momentum along the proton beam direction), arising in central exclusive diffractive processes 
[13, 3], provides a unique possibility to study directly the coupling of the Higgs-like bosons to 
the bottom quarks, because the LO QCD background is strongly suppressed by the Jz = rule. 
As is well known, the determination of the Hbb Yukawa coupling appears to be very difficult 
for other search channels at the LHC. 

The theoretical formalism [14, 1, 15] needed to describe a central exclusive diffractive process 
of a system A contains quite distinct parts, shown symbolically in Fig. 1. In brief, we first have 
to calculate the gg ^ A hard subprocess, H, convoluted with the gluon distributions fg. Next 
we must account for the higher loop corrections which refiect the absence of additional QCD 
radiation in the hard subprocess - that is, for the Sudakov suppression factor T. Finally we 
must enter soft physics to calculate the survival probability S"^ of the rapidity gaps either side 
of ^ - that is the probability that the exclusive nature of the process will not be destroyed by 
the secondaries produced by the rescattering of the incoming particles. 

2 Objective 

The uncertainties associated with the prediction of the rate of a central exclusive process are 
potentially not small. Each of the above stages has its own uncertainties. Therefore, it is 
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important to perform checks of the approach using processes with appreciable cross sections 
that will be experimentally accessible in the first data nms of the LHC with integrated lu- 
minosities in the range 100 pb"^ to 1 fh-\ Here, our aim is to identify processes where the 
different ingredients of the formalism used to calculate central exclusive production can be 
tested experimentally, more or less independently. 

The outgoing protons in the forward regions of the main LHC detectors will be measured 
by the proton tagging detectors (Roman Pots). At 220 m on either side of the CMS detector 
there exists the roman pots of the TOTEM experiment [16, 17]. ATLAS will have the ALFA 
detector [18] at 240m as well as the proposed RP220 detector [19] at 220m. As we already 
mentioned, FP420 project [8, 6] proposes to install forward proton taggers for both ATLAS 
and CMS. Note that RP220 and FP420 aim to operate at high luminosities. However it is 
quite likely that for the first 1-2 years of LHC running the forward proton detectors will not 
be operational. We therefore first consider measurements which do not rely on tagging the 
forward outgoing protons. 

Even without proton tagging, diffractive measurements at the LHC can be performed 
through the detection of rapidity gaps. This is a well known technique used extensively at 
HERA and the Tcvatron. A summary of the forward detectors instrumented around ATLAS 
and CMS is given, for instance, in [20]. The central detectors (CD) of CMS and ATLAS have 
an acceptance in pseudorapidity rj of roughly Ir^l < 2.5 for tracking information and |?7| < 5 for 
calorimeter information. We will discuss a situation where a heavy system A is detected in the 
region \r]\ < 2.5, and where the calorimeters in 2 < 77 < 5 interval, are used to select events 
with rapidity gaps. 

In the present paper the word gap means a rapidity interval Arj devoid of hadronic activity 
- for the charged particles we assume a 'track veto', while the absence of neutrals should be 
checked by the calorimeters. The amplitudes of all the processes discussed below are infrared 
stable and are not affected by the possibility of soft gluon emission. 

Due to angular ordering, which originates from coherence, in those processes where the 
gap is provided by 1^-boson exchange (as in Fig. 2(a,c) below), the only possible soft gluons 
occur at the edge of the gap and arise from the corresponding quark jet. These gluons should 
be accounted for in the jet searching algorithm. In cases where the gaps are associated with 
colour-singlet, two-gluon exchange (as in central exclusive dijet production), the presence of 
Sudakov-like T-factors in the unintegrated gluon distributions (see eq.(7) below) guarantee that 
there is no emission of any additional soft gluons. 

The selection of rapidity gap events by a 'veto' trigger can be used up to rather large 
luminosities, when the mean number, A^, of interactions per bunch crossing is sizeable. However 
at larger luminosities the efficiency of the trigger is reduced by a factor ~ that is, by the 
probability to have no additional 'pile-up' inelastic interaction in the bunch crossing. This 
probability can be measured independently in the same experiment. 

Of course, the proposal to use calorimeters in the 2 < rj < 5 interval to select the events with 
a rapidity gap does not mean that we will only consider gaps with Ar] < 3. First, part of the gap 
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can be at smaller 77 and, secondly, extensive additions are foreseen to enlarge the coverage in the 
forward regions^. Possible "holes" in particle observation in small rapidity intervals between 
different calorimeters do not affect our predictions very much. The probability to produce 
extra soft hadrons in the hole region only is proportional to the effective Pomeron-Pomeron 
cross section which is rather small, according to the triple- Regge analysis of UA8 data [22]. We 
evaluate the correction to be less than 3% for a hole of size A^^hoic ~ 1-5. 

The main uncertainties of the predictions for exclusive processes are associated with the 
calculation of 

(i) the probability S'^ that additional soft secondaries will not populate the gaps separating 
the centrally produced system A from the outgoing protons (or the products of their 
dissociations)^; 

(ii) the probability to find the appropriate gluons that are given by generalized unintegrated 
distributions fg{x, x',Qf); 

(iii) the higher order QCD corrections to the hard subprocess amplitude, where the most 
important is the so-called Sudakov suppression caused (in the Feynman gauge) by the 
double-log loop, denoted T in Fig. 1, see [14]; 

(iv) the so-called semi-enhanced absorptive corrections (see, for example, [26, 27, 28]) and 
other effects which may violate the soft-hard factorization (which is imphed by Fig. 1). 

We discuss below, in turn, possible checks of the various ingredients of the calculation of 
the exclusive cross sections. To be precise, we address the uncertainties (i)-(iv) in Sections 3-6 
respectively. 

^CMS will have the CASTOR calorimeter operating with the coverage of 5.1 < \r]\ < 6.5 and a Zero 
Degree Calorimeter (ZDC) with an acceptance for neutral particles with ry > 8. Both calorimeters have an 
electromagnetic and hadronic section. Moreover, CMS is expected to have integrated read out with TOTEM 
(sec [16]) allowing CMS to benefit from the TOTEM forward coverage and TOTEM, in turn, from the CMS 
central coverage [17]. ATLAS plan a Cerenkov detector LUCID with an acceptance 5.4 < 77 < 6.1 and a ZDC 
with an acceptance 8.3 <r] < 9.2. Diffractive studies are under discussion also at ALICE, see for instance, [21]. 
The ALICE detector has a central barrel covering the pseudorapidity range —0.9 < r] < 0.9 and (on one side) 
a muon spectrometer covering the region of 2.4 < r) < 4 and a ZDC. Additional detectors for trigger purposes 
and for event classification are placed on both sides of the central barrel, such that the range — 3.7<77<5is 
covered. This configuration allows the possibility of a (double) rapidity gap trigger by requiring no activity in 
the event classification detectors, sec [21]. 

^There are quite a few theoretical studies of the S"^ factor, starting from the publications [23, 24] to more 
recent ones [15, 25] and references therein. Following Bjorken ([24]) this quantity is often called the survival 
probability or soft survival factor. 
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Figure 2: Diagrams for (a) W production with 2 rapidity gaps, (b) inclusive W production, 
and (c) Z production with 2 rapidity gaps. 



3 Gap survival probability 5*^ 

As a rule, the gap survival probability, that is the probability that the secondaries produced in 
additional soft interactions do not populate the rapidity gaps, is calculated within a multichan- 
nel eikonal approach. This method of calculation implicitly assumes a factorization between 
the soft and hard parts of the process. The probability of elastic pp rescattering, shown sym- 
bolically by S"^ in Fig. 1, can be calculated in a model independent way once the elastic pp 
amplitude is known, that is, the elastic cross section daei/dt is measured at the LHC. However 
there may be some excited states (corresponding to A^* resonances and low mass diffractive 
dissociation) between the blob S and the perturbative QCD amplitude on the right-hand-side 
of Fig. 1. The presence of such additional states enlarges the absorptive correction. Usually this 
effect is calculated using a Good- Walker formalism [29] with a two- or three-channel eikonal. In 
order to experimentally check the role of this effect, we need to consider a process with a bare 
cross section that can be reliably calculated theoretically. Good candidates are the production 
ofWorZ bosons with rapidity gaps on either side. 



3.1 W production with rapidity gaps 

In the case of W+gaps' production the main contribution comes from the diagram shown in 
Fig. 2(a) [30]. One gap, A?7i, is associated with photon exchange, while the other, , is 
associated with W exchange. The cross section of this process can be calculated straightfor- 
wardly. It is proportional to the quark distribution in the proton at a large scale and not too 
small X, where the uncertainties of the parton densities are small. To select these events at the 
expected luminosities of the early data runs, we can use the rapidity gap veto trigger combined 
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with a high pt decay lepton or jet trigger. The momentum transferred through the photon is 
typically small and the rapidity gap A771 runs from the rapidity of the W boson, yw, up to the 
maximum rapidity measured. The gap, corresponding to W exchange, is limited by the 
rapidity of the quark jet, which balances the transverse momentum of the W boson (which has 
a broad distribution with qtw ~ Mvk/2 since it is driven by the ^-channel W propagator). 

At first sight, the probability for soft rcscattering in such a process is rather small and we 
would expect S*^ ~ 1. The reason is that the transverse momentum, kt, distribution of the 
exchanged photon is given by the logarithmic integral 



for which the dominant contribution comes from the low /c| region. In other words, the main 
contribution comes from the region of large impact parameters, bt, where the opacity of the 
proton is small. However the minimum value of \t\, 



is not negligibly small. Note that the momentum fraction Xp = l — ^ associated with the upper 
proton can be measured with sufficient accuracy^, even without the tagging of the forward 
protons, by summing the momentum fractions 



of the outgoing W and the hadrons observed in the calorimeters. As long as the gap Ari2 
is large, the dominant contribution to the sum C — comes from the decay products of 
the W boson^. For example, for rjw — 2.3(— 2.3), we expect an ^ distribution centred about 



The predictions are shown in Figs. 3-5. First, in Fig. 3, we neglect the absorptive correction, 
and show the cross section for W+gaps events for different choices of the gap size A772. We see 
cross sections of up to 1 pb, depending on the gap size and on The rcscattering decreases 
the cross section by the factor S'^ shown in Figs. 4 and 5. Already, at this stage, we face the 
possibility of a violation of soft-hard factorization. We use the Good- Walker formalism [29] to 
describe multichannel eikonal rcscattering. That is, we use the diffractive eigenstates (f)k which 
diagonalize the 'nucleon'-Pomeron couplings jSij, which describe the transition from nucleon 
excited state i to state j. In Ref. [15] it was assumed that each eigenstate has the same size 
and the same partonic composition. The continuous curve in Fig. 4 was calculated under this 
assumption. On the other hand, in the extreme cases where we assume that the partons which 
participate in the process of Fig. 2(a) are concentrated in the diffractive eigenstate with the 

^The CDF collaboration [31] have demonstrated that this method provides an accurate determination of ^. 
^Evcn for th(^ Icptonic decay channels, the momentum of the W boson can be reconstructed with knowledge 
of the missing transverse energy and Mw- 
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C ~ 0.1(0.001). 
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Figure 3: The cross sections for ly+gaps events as a function of ^ at the LHC for different 
choices of the size of the rapidity gap A?72- No suppression from rescattering effects is included. 
The momentum fraction of the proton that is carried by the photon is ^ = 1 — , where Xp is 
the momentum fraction of the outgoing upper proton that emits the photon in Fig. 2(a). We 
use the superscript + to indicate that the cross sections correspond to the configuration where 
the photon is emitted by the proton going in the + direction. The event rate observed in the 
central detector will be twice as large as that shown. 
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Figure 4: The rapidity gap survival factor 5*^ as a function of C, calculated using the global soft 
model of [15]. 



lowest (largest) absorptive cross section, we obtain the results shown by the dashed (dotted) 
curves. A rather more realistic scenario (see, for example, [26]) is to allocate the valence 
quark to the component with the smallest absorption and the sea quark to the component with 
largest cross section. The corresponding predictions are shown in Fig. 5. Since the valence 
quark contribution is more important for production and for the configuration with the 
largest gap size the expected gap survival factor S'^ is found to be larger. Now, the spread 
of predictions caused by the different partonic content, shown in Fig. 5, is rather small. With 
a more realistic assumption (in which each diffractive state contains some part of the valence 
and some part of the sea) the spread will be even smaller. 

In the first LHC data runs it may be difficult to measure the absolute value of the cross sec- 
tion with sufficient accuracy. Most probably the ratio (W+gaps/W inclusive) will be measured 
first. In this case, the inclusive W production process (in the same kinematic region. Fig. 2(b)) 
plays the role of the luminosity monitor. Note that the cross section for inclusive W production 
is much larger than that with rapidity gaps. The reason is an inclusive W is produced directly 
by qq fusion, which is prohibited for gap events since the colour flow produced by the t-channel 
quarks populates the 2 rapidity gaps. 

Of course, the survival factor S"^ measured in ly+gaps events is quite different from that 
for exclusive Higgs production, which comes from smaller values of bt. Nevertheless this mea- 
surement is a useful check of the model for soft rescattering. 
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Figure 5: The rapidity gap survival factor 5^ as a function of ^ calculated using the global soft 
model of [15], assuming that the valence (sea) quarks are associated with the weak (strong) 
absorptive components. The small spread of the predictions arising from the different partonic 
content of the diffractive eigenstates mean that 1^+gaps events offer a meaningful test of the S"^ 
factor. Note that S"^ for the W'^ signal is larger since it has a bigger valence quark contribution. 
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3.2 Z production with rapidity gaps 

It would appear that a good way to study the low region directly is to observe Z boson 
production via WW fusion, see Fig. 2(c). Here, both of the rapidity gaps originate from heavy 
boson exchange and the corresponding ht region is similar to that for central exclusive Higgs 
production^. The expected Z+gaps cross section is of the order of 0.2 pb, and 5'^=0.3 for 
A77i,2 > 3 and for quark jets with Et > 50 GeV [32] . Again it would be sufficient to measure 
the ratio {Z+gaps/Z inclusive). 

One problem is that even with the Et > 50 GeV cut, the QCD background arising from the 
QCD bb central exclusive production is comparable to the electroweak — > Z + 2 jet signal. 
Therefore we should concentrate on the leptonic decay modes of the Z boson, which results in 
a smaller event rate^. 



3.3 With tagged protons 

So far we have discussed measurements which do not depend on observing the forward protons. 
However, when the forward proton detectors become operational we can do more. Both the 
longitudinal and transverse momentum of the forward protons can, in principle, be measured. 
Hence we can study the kt behaviour of the cross section section for 1^+gap events, and scan 
the proton opacity, as described in [30]. 



4 Generalized, unintegrated gluon distribution fg 

The cross section for the central exclusive production of a system A is calculated using Qt 
factorization, and essentially has the form [14] 



^2 

a{pp^p + A + p) ~ — 



SJ Q 



a{gg ^ A). (4) 



The first factor, 5"^, is the soft survival factor discussed in the previous section and the factor 
1/6^ arises from the integration over the transverse momentum of the forward proton assuming 



^Even so it may not be exactly the same bt region, since, in general, the impact parameter distributions of 
quarks and gluons may be different. 

^Note that in the recent study [33] it was demonstrated that the so-called Track Counting Veto (TCV) is 
robust for selection of the central rapidity gap events in vector-boson fusion H t^t" searches at CMS. 
The idea of track counting is close in spirit to the hadron level gap selection of [23, 32] and, as compared to 
the (more) standard calorimeter jet veto technique, the TCV has various advantages. In particular, it does 
not involve calorimeter scale uncertainties, suffers less from pile-up contributions and so can be used at higher 
luminosities. Moreover, it is planned to apply this technique at CMS using data on central Z IJi'^ H~ events. 
Though the original motivation for this analysis [33] is to study the detector conditions, such a measurement, 
in particular, presented in terms of the ratio (Z+gaps/Z inclusive) may appear to be one of the first tests of 
models for soft survival factor. 
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the form da /dpi oc exp(— Also, fg denotes the generahzed, unintegrated gluon distribution 
in the hmit of pt — > 0. Below we give a more precise definition of fg and explain how it is 
determined. 

The generalized gluon distribution has not yet been measured explicitly. However, for the 
case of interest, where the skewedness is small, it can be obtained from the conventional diagonal 
gluon distribution, known from the global parton analyses. The procedure is as follows. We 
consider the skewedness first in the transverse momenta, and then in the longitudinal momenta, 
carried by the two t-channel gluons in Fig. 1. The transverse momentum, pt, transferred through 
the rapidity gap in a central exclusive process is limited by the proton form factor. This is not 
a large transverse momentum: p| ~ — t <S Qj. Typically, we have the hierarchy 

l^minl < 10-" GeV^; p2_0.2GeV2; _ 4 GeV^; (5) 

which follows, respectively, from (2) with the relevant ^ ~ 0.01 or less; from the proton form 
factor; and, finally, from the presence of the Sudakov-like form factor, T of (8), in the integrand 
of (4) - the maximum of the integrand occurs at ~ 4 GeV^, see [1, 13]. 

As a consequence, it is reasonable to assume the following factorization of the gluon dis- 
tribution: F(t) fg{x,x' ,Ql, fi'^). That is, we work in terms of the distribution fg{x,x' ,Ql, fi'^), 
which is skewed only in longitudinal momentum fractions, a; ^ a;' in Fig. 1; is essentially 
the transverse momentum of each i-channel gluon. Finally, /x^ ~ is the factorization 

scale which separates the partonic distribution from the matrix element of the hard subprocess 
gg A. The precise scale, to be used, is given in (13) below. 

In the region of interest 

[X —=\ < X ~ < 1. (6) 

The generalized distribution can be obtained using the Shuvaev transform [34] which assumes 
that the Gegenbauer moments of the generalized parton distributions are equal to the Mellin 
moments of the diagonal distributions. The accuracy of this assumption is 0{x'^). As a result 
we obtain [35]^. 

fg{x, x', Ql, I?) = Rq q^^ q2 {^T{Qt,ii) xg{x, Q^)^ , (7) 

where Rg accounts for the x ^ x' skewedness. Note that the double log Sudakov suppression 
T{Qt,fi) is now included in the unintegrated gluon distribution fg, since to provide Qt fac- 
torization we choose an axial/planar gauge (and not the Feynman gauge used to draw T in 
Fig. 1). The Sudakov form factor T ensures that the gluon remains untouched in the evolution 
up to the hard scale so that the rapidity gaps survive. It results from resumming the virtual 
contributions in the DGLAP evolution, and is given by [36] 



T{Qt, 1^) ^ exp ' ' ^^t' ■ 



27r k 



2 



dz], (8) 



''In the actual computations we use a more precise form as given by eq.(26) of Ref.[35] 
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(a) (b) 




Figure 6: Exclusive T production via (a) photon exchange, and (b) via odderon exchange. 



where A is specified in (13) below. 

The main uncertainty here comes from the lack of knowledge of the integrated gluon dis- 
tribution g{x, Qf) at low x and small scales. For example, taking Qf = A GeV^ we find that a 
variety of recent MRST [37] and CTEQ [38] global analyses give gluons which have a spread of 

xg = {3- 3.8) for X = 10"^ and xg = (3.4 - 4.5) for x = 10^1 (9) 

These are big uncertainties bearing in mind that the cross section for central exclusive produc- 
tion depends on (xg)^. 

4.1 Exclusive T production as a probe of fg 

To reduce the uncertainty associated with fg we can measure exclusive T production. The 
process is shown in Fig. 6(a). The cross section for jp —>■ Tp [39] is given in terms of exactly 
the same generalized, unintegrated gluon distribution fg that occurs in Fig. 1. If the T is 
detected at the LHC in the rapidity interval {\y\ < 2.5) then it will sample the Xi^2 intervals 
(10~^ — 10^^). To probe even lower x we can either detect more forward T production at 
ALICE or LHCb [40] or study exclusive J/ip production, see [30]. 

Of course, there may be competition between production via photon exchange. Fig. 6(a), 
and via odderon exchange, where the upper three t-channel gluons form the C = — 1 odderon 
state^, see Fig. 6(b). To date, odderon exchange has not been observed. On the other hand, 
a lowest order perturbative QCD calculation indicates that the odderon process (b) may be 
comparable to the photon- initiated process (a), see, for example, [41]. If the upper proton is 
tagged, it will be straightforward to separate the two mechanisms since odderon production has 
no 1/g^ singularity characteristic of the photon. Photon exchange populates mainly the low qt 
region, that is at low pt of the outgoing proton, whereas odderon production occurs at relatively 

®Here the colour indices of the gluons are convoluted with the dabc tensor. 
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large pt- Without a proton tagger, an odderon would be revealed by a measured cross section 
that is larger than that predicted for photon exchange. Exclusive T (or J/ip) production could 
be the first hint of the odderon's existence. 



Returning to process (a), we note that the cross section is of the form 



a 



(a) 



(10) 



where the photon flux Nj is well known. For small qt, we may neglect the proton form factor 
and use the leading log approximation 



at the LHC energy^. The signal will be diluted by the T ^ /i/i branching fraction of 0.025. 

In order to use this process to constrain the gluon distribution of the proton it would 
be preferable to tag the lower proton. Otherwise there will be some admixture of proton 
excitations and to calculate precisely the gap survival factor S'^, as before, we will have to 
resolve the partonic content of the different diffractivc eigenstatcs. On the other hand, for low 
values of ^, the expected value S'^ is close to 1 (sec the continuous curve in Fig. 4), and it may 
be sufficient to use the existing HERA data for the ratio of the cross sections of diffractive J/ t/j 
photoproduction with and without the proton dissociation, to extract the contribution to the 
cross section from 'elastic' events in which the lower proton does not dissociate. 

5 Improvements to the hard subprocess 

The hard matrix elements necessary to predict the cross section for a central exclusive process 
are calculated using perturbative QCD mainly at leading log accuracy [1]. Of course, there 
exist numerous NLO results, but as a rule these are applicable to inclusive processes which do 
not allow for the fact that the centrally produced system A must be in a colour singlet state 
and obey the Jz — selection rule, that is the projection of the total angular momentum along 
the incoming proton directions should be zero. Moreover the relation between the generalized, 
unintegrated gluon distribution, fg, and the conventional diagonal gluon, g, was also based on 
a LO calculation [35]. 

^The estimate assumes that the forward protons are not tagged, and so is enhanced by a factor of about 2 
caused by the dissociation of the lower proton in Fig. 6(a). 




(11) 



dy 



y=o 



~ 50 pb 



(12) 
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Numerically, the largest effects may come from the next-to-leading-log (NLL) corrections to 
the double log Sudakov-likc factor T of (8). To include these NLL corrections phenomenolog- 
ically we choose the limits of the integration over /c^ in such a way to reproduce the result of 
an explicit first- loop calculation. It has been shown (see footnote 9 in [9]) that this is achieved 
by choosing 

k 

A = V with u = 0.62Ma (13) 

pt + h 

in (8). 

On the other hand, the lower limit of the integration has not yet been validated at 
this level. Clearly the contribution from the very small kt region vanishes due to destructive 
interference between the emissions from the active gluon (x) and the screening gluon {x'). 
However, variation of the lower limit from Qt/2 to 2Qf may alter the prediction by up to an 
order of magnitude. 

Fortunately, the contribution from the region of kt ~ Qt can be calculated with better 
precision. We are seeking the next-to-leading log correction to the double logarithmic T factor 
of (8), that is, for the single log terms in T. Therefore, in the region kt ~ Qt we need to 
keep only the LO BFKL-like term which contains the longitudinal logs. This term is given by 
the usual BFKL kernel [42] which sums all the leading a^logMA contributions. If we consider 
the kt integration in (8) with a lower cut-off of k^, then the above summation amounts to the 
replacement 

ky ^ /c| 1^ kl + {Qt-kt)V"' JQl kl"" 

First, we notice in the middle expression for the integral above that for kt ^ Qt the last term 
in the brackets is negligible, while for kt ^ Qt the whole expression in brackets goes to zero 
and removes the infrared divergence. Now we can put the infrared cut-off ko = 0. The result 
of an explicit calculation of the middle expression for the integral shows that it is equivalent to 
the integral in (8) with the lower limit given by Qt, as indicated in (14). 

Thus, in summary, both the upper and lower limits of the kf integral in (8) are fixed so as 
to reproduce the one-loop contributions. The potentially large ambiguity in the value of T has 
been removed. 

To check experimentally the formalism used in the perturbative QCD calculations for the 
central exclusive matrix element we should study an exclusive process with the emission of one 
additional jet. Note that the physical origin of the Sudakov T factor is that it expresses the 
probability not to emit additional gluons. The formula for T{Qt,iJ.) of (8) can be written as 
T — exp(— n) where n is the mean number of gluons emitted in the interval {QtiH). Thus, 
the observation of the explicit emission of additional gluon jets would provide a direct test of 
the formalism. Since the central system must be colour neutral, we need either consider the 
emission of two extra jets or to have the possibility to compensate the colour of one emitted 
gluon by rearrangement of the colour content of the system A. The optimum choice is to 
observe the emission of a third jet in the production of a pair of high Et jets, that is, in high 
Et dijet production. 
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Figure 7: The rapidities of the three jets in the central system. Note that the rapidity ua of 
the whole central system does not necessarily occur at ?/ = 0. The rapidity interval containing 
the three jets is denoted by 5rj, outside of which there is no hadronic activity. 



5.1 Three-jet events as a probe of the Sudakov factor 

Traditionally, the search for the exclusive dijet signal at the Tevatron, pp p + jj + p, is 
performed [43] by plotting the cross section in terms of the variable 

Rn = M,,/Ma . (15) 

Ideally, we might expect exclusive dijet production to show up as a narrow peak centred at 
Rjj = 1, since, for these events, the mass of the dijet system, Mjj, is equal to the mass. Ma, of 
the whole central system. Unfortunately, in practice, the Rjj distribution is strongly smeared 
out by QCD bremsstrahlung, hadronization, the jet searching algorithm and other experimental 
effects [43, 44]. 

To weaken the role of this smearing it was proposed in Ref. [44] that the observed dijet 
distribution be studied in terms of a new variable 

Rj = 2Et (cosh ri*)/MA , (16) 

where only the transverse energy Ex and the rapidity r] of the jet with the largest Et are used 
in the numerator. Here Tj* = rj — ua where ua is the rapidity of the whole central system^". 
Clearly the jet with the largest Et is less affected by hadronization, final parton radiation etc. 
In particular, final state radiation at lowest order in as will not affect Rj at all, since it does 

^"Note that the transverse momentum of the dijet system can be neglected, since it is very small compared 
to the Et resolution. 
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not change the kinematics of the highest Et jet used to evaluate (16). Even with the emission 
of an extra jet during the final parton shower, we will still have Rj — 1. Thus, to see the role of 
QCD radiation on the Rj distribution, we only have to account explicitly for additional gluon 
radiation in the initial state. At leading order, it is sufficient to consider the emission of a third 
gluon jet, as shown in Fig. 7, where we take all three jets to lie in a specified rapidity interval 
Srj. The reason why it is sufficient to consider only one extra jet, is that the effect of the other 
jets, which, at LO, carry lower energy due to the strong ordering, is almost negligible in terms 
of the Rj distribution. 

The cross section da/dRj, as a function of Rj, for the exclusive production of a high Et 
dijet system accompanied by a third (lower Et) jet was calculated and discussed in detail in 
[44]. The mass of the central system, M^, can be accurately measured by the missing mass 
reconstructed from the tagged protons, or, failing that, by accurately summing both the light- 
cone momentum fractions ^'^ and of the hadrons observed in the calorimeter, see (3). As 
shown in [11], the usage of the Rj variable is quite beneffcial for the extraction of exclusive dijet 
events. In particular, the dependence on the jet selection criteria (for example, the cone radius 
and it!-parameter) is less marked for the Rj than for the conventional Rjj variable. Further 
discussion can be found in [44]. Moreover, note that studying the Rj distribution, we do not 
need to select events with a rather large Et of the third jet. To calculate Rj it is sufficient to 
measure only the Et of the largest -Et jet, together with the rapidities of first and second jets, 
and the values of Ma and ua as measured via ^'^ and 

In Fig. 8 we show the Rj distribution of both exclusive two- and three-jet production 
expected at the LHC. For three-jet production we show predictions for two choices of the 
rapidity interval 57] within which the third jet must lie. If we take the largest Et jet to have 
Et > 50 GeV at the LHC, we see that the cross section for exclusive three-jet production reaches 
a value of the order of 100 pb. Of course, if we enlarge the rapidity interval 5r] where we allow 
emission of the third jet, then da/dRj will increase, see Fig. 8. Indeed, the measurement of the 
exclusive two- and three-jet cross sections as a function of Et of the highest jet allow a check 
of the Sudakov factor; with much more information coming from the observation of the 5r] 
dependence of three-jet production. Note that the background from double-Pomcron-cxchange 
should be small for Rj > 0.5, and can be removed entirely by imposing an Et cut on the third 
jet, say Et > 5 GeV. 

Another way to observe the effect of the Sudakov suppression is just to study the Et 
dependence of exclusive dijet production. On dimensional grounds we would expect da/dE^ oc 
1/E^. This behaviour is modified by the anomalous dimension of the gluon and by a stronger 
Sudakov suppression with increasing Et- Already the existing CDF exclusive dijet data [43] 
exclude predictions which omit the Sudakov effect. In the Et interval from 10 to 35 GeV the 
Tevatron observations show that the cross section falls by an order of magnitude faster than 
the prediction [45] , based on [46] , which does not include the Sudakov suppression, see Fig. 20a 
of [43]. 

It is clear that precise measurements of the Et behaviour of the exclusive dijet cross section 
at the LHC offer the possibility to probe the corrections to (8). However the study of exclusive 
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Figure 8: The Rj distribution of exclusive two- and three-jet production at the LHC. Without 
smearing, exclusive two-jet production would be just a (5-function at Rj = 1. The distribution 
for three-jet production is shown for two choices of the rapidity interval, 6t], containing the jets; 
these distributions are shown with and without smearing. Here, we have taken the highest 
jet to have Et > 50 GeV. To indicate the effect of jet smearing, we have assumed a Gaussian 
distribution with a typical resolution a = O.G/y/Er in GeV. 
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(a) (b) (c) 




Figure 9: (a) A typical enhanced diagram, where the shaded boxes symbohcally denote fg, and 
the soft rescattering is on an intermediate parton, giving rise to a gap survival factor ^cn; (b) 
and (c) are the Reggeon and QCD representations, respectively. 

three-jet events will provide much more information and, indeed, allow a study of the integrand 
of (8). 

6 Soft-hard factorization: enhanced absorptive effects 

We mentioned in Section 3 that the soft-hard factorization implied by Fig. 1 may be already 
violated if the different diffractive eigenstates of the pp interaction have different partonic 
distributions. To date, there is no unambiguous model to distribute the partons obtained in 
the global analyses between the different diffractive components. However, this is not a major 
uncertainty. It is seen in Fig. 5 that the difference between and W~, which have a different 
valence quark contribution, is not large. 

Another potential source of the violation of soft-hard factorization arises from the so-called 
enhanced Reggeon diagrams, which occur from the rescattering of an intermediate parton gen- 
erated in the evolution of fg. Such a diagram is shown in Fig. 9(a). Fig. 9(b) is the Pomeron 
representation of the diagram, where the rapidity of the intermediate parton fixes the position 
of the triple-Pomeron vertex. Clearly such diagrams will violate the soft-hard factorization. 

The contribution of the first Pomeron loop diagram. Fig. 9(b) was calculated in perturbative 
QCD in Ref. [27]. A typical perturbative diagram is shown in Fig. 9(c). For LHC energies it 
was found that the probability of such rescattering may be numerically large^^. The reason 

^^Of course, the higher-order Reggeon diagrams will considerably reduce the size of the effect. Moreover, 
allowing for the more complicated multi-Pomeron vertices makes the suppression of the first Pomeron loop 
contribution even stronger. This effect has be seen in a particular model, see Fig. 17 of [47]. It is practically 
impossible to quantify the reduction within the framework of perturbative QCD. 
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is that the gluon density grows in the low x region and, for low kt partons, approaches the 
saturation limit. The recent estimates of [48] indicate that in the region of x~ ~ 10~^ (which is 
relevant for the LHC kinematics) the value of the saturation scale, Qs, exceeds 1 GeV. In other 
words, a parton with kt < Qs will receive an important absorptive correction. However, this is 
true for the central region of impact parameter space. On the other hand, eikonal absorption, 
which is shown in Fig. 1 and discussed in Section 3, makes the centre of the disk almost black. 
This eikonal factor, S"^, already forces the central exclusive signal to occur at relatively large 
bt, namely bt > 0.5 — 0.6 fm, see Fig. 22 of [47]. In this peripheral region of the proton, the 
value of Qs is rather low. In fact, (5| ~ 0.3 GeV^ is found in Fig. 11 of [48]. As a consequence, 
the enhanced diagram will affect only the very beginning of the QCD evolution - the region 
that cannot be described perturbatively and which, in our calculation of the cross sections for 
central exclusive processes, is already included phenomenologically - that is, described in terms 
of the multichannel eikonal framework. 

This is probably one reason why semi-enhanced screening corrections are not seen in leading 
neutron production measured at HERA [49] . As the energy of the virtual photon increases, the 
available rapidity region grows, so the number of intermediate partons that may participate in 
the rescattering increases, and we expect that the fraction of leading neutrons should decrease. 
This is not seen in the data [49]. Further discussion is given in [26, 28, 50]. 



6.1 Pomeron-p interactions as a probe of the enhanced effect 

Experimentally, it appears at first sight, that we may study the role of semi-enhanced absorp- 
tion by observing the l^-|-2 gaps process shown in Fig. 2(a). To do this one may vary the 
transverse momentum of the accompanying quark jet, q, and the size of the rapidity gap A772. 
For lower transverse momentum of the quark jet we expect a stronger absorptive effect, that 
should decrease with increasing A?72, since the number of partons in the rest of the rapidity in- 
terval increases. Unfortunately, the photon-initiated process of Fig. 2(a) occurs at large impact 
parameter bt where the probability of rescattering is small. Moreover it will be very difficult 
to observe the quark jet at low Et in the main calorimeter (—5 < rj < —3). Therefore we will 
discuss other processes which depend on rescattering on intermediate partons. 

The observations we have in mind are the measurements of ratio R of diffractive (one-gap) 
events for W (or T or dijet) production as compared to the number of events for the inclusive 
process (shown in Fig. 2(b) for W production). These processes are shown schematically in 
Fig. 10. In other words, R is the ratio of the process in diagram (c) to that in diagram (a). 
That is 

^ no. of {A + gap) events ^ a'^'^{x]p, (3, fi^) /c,2c2v /.yx 
no. of (inclusive A) events c^-'^^x = Pxjp,fi^) ^ ^"^^""'"^ ^ ' 

where a"^'^' and a'^^^ are the parton densities determined from the global analyses of inclusive 
and diffractive deep inelastic scattering data, respectively. The heavy central system A is either 



18 



(a) (b) (c) 




Figure 10: Schematic diagrams for (a) the inclusive production of a system A, (b) and (c) for 
the diffractive production of A without and with 'enhanced' soft rescattering on intermediate 
partons. The system A is taken to be either a W boson or an T or a pair of high Et jets. 

W or a pair of high Et jets or T or a Drell-Yan yU+yU" pair. For W or yU+yU" pair production 
the parton densities a are quark distributions, whereas for dijet or T production they are 
mainly gluon densities. The diffractive parton densities are known from analyses [51, 52] of 
diffractive deep inelastic data. Thus measurements of the ratio R will probe the gap survival 
factor averaged over the impact parameter bt- 

If we neglect the effect of enhanced rescattering on intermediate partons, that is set S"^^ = 1 
then R is the ratio of process (b) to (a). However this ratio will be reduced by the rescattering 
on the intermediate partons as sketched in Fig. 10(c)). The largest contribution to the extra 
suppression factor S'^^ comes from the interaction of the upper intermediate partons with the 
lower proton, since the partial energy corresponding to this interaction is larger, which leads 
to a larger absorptive cross section. Note that the enhanced effect is practically forbidden in 
a rapidity interval close to the rapidity of A (shown by the bold lines in Fig. 10(c)), as the 
transverse momenta of partons in this region are quite large. The mean number of partons, 
(A^), emitted in the evolution from 1 GeV up to the hard scale yU ~ Ma/2 is given by the 
power of the exponent in the T factor of (8). We have (N) ~ 1 or 2 for Ma — 10 or 100 GeV 
respectively, that is for exclusive T or ly production. These are the partons with pt > Qs{bt) 
which do not suffer the enhanced absorptive effect. As a rule, each parton occupies about 
one unit of rapidity. Therefore in the evaluation of the enhanced effect we exclude from our 
calculation an appropriate rapidity interval on either side of A, which is shown by the bold lines 
in Fig. 10(c). The presence of a "threshold" factor, that is the existence of a rapidity interval 
on either side of A where the enhanced contribution is effectively forbidden, was emphasized in 
[26]. An analogous threshold strongly suppresses the enhanced absorption to exclusive Higgs 
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production [28]. 

In the region of relatively large bt the value of Qs is low. Hence we evaluate the enhanced 
absorption in terms of Reggeon framework, using the phenomenological triple- Pomeron coupling 
extracted from diffractive J/ip production [53]. The absorptive corrections for J/ip production 
are low, and the value g^p ~ Qn/^ should be close to the original bare triple- Pomeron coupling. 
Here, is the coupling of the Pomeron to the proton. The contribution of the one-loop Regge 
diagram of Fig. 9(b) can be quite large. Indeed, its strength, as compared to the original cross 
section, can be estimated by evaluating Fig. 9(b). We obtain 



where crpp{r)') is the proton-proton cross section evaluated at the energy y/s — mjv exp(r7'/2). 
The factor 4 comes from the AGK cutting rules [54] and S ~ 5 GeV~^ is the sum of the t slopes 
of the form factors involved. The integral is taken over the appropriate rapidity intervals of the 
intermediate partons in Fig. 10(c). Since the result is sizeable, we allow for the summation of 
the higher order diagrams by replacing the absorptive factor 1 — r by 5*^^ = cxp(— r). 

Experimentally, we can measure the rapidity ua of the central system {A = W, T or dijet) 
and also the momentum fraction carried in the gap direction by the soft hadrons, ^~ = J2d, 
recall the analogous sum of (3). Thus we know the value of xjp = ^~ + ^J- That is, we can 
observe a double distribution (fa'^^^/dxipdyA, and form the ratio R using the inclusive cross 
section, da^^^^/dyA- If we neglect the enhanced absorption, it is straightforward to calculate the 
ratio R of (17). The results are shown by the dashed curves in Figs. 11 and 12 as a function 
of the rapidity yA of the heavy system A. When we allow for the enhanced rescattering the 
ratios are reduced, and lead to steeper yA distributions, as shown by the continuous curves. 
These plots correspond to a fixed value of ^~ = 10 GeV/7 TcV. That is, we assume that the 
soft hadrons observed in the central calorimeter carry a longitudinal momentum^^ p^^lO GeV 
in the direction of the rapidity gap associated with Pomeron exchange^^. 

The results shown in Figs. 11 and 12 should be regarded as an indication of the size of 
possible enhanced effects and not as quantitative predictions. First, the model used to estimate 
5*^^ is quite naive. Next, the results for large negative yA sample the diffractive gluon density 
at large /?, where it is not well constrained. Finally, for T (and dijet production at the lower 
Et values) we need, for large positive yA, the conventional gluon density at very small x — 
(3x]p ~ M^exp(— t/^)/i/i ~ 10~^, where it is not well determined in the global analyses. 

First, we discuss the results shown in Fig. 11 for the ratio R for W production at positive yA- 
Here, due to the large W mass, we have no available rapidity interval for enhanced rescattering, 

^^This value oipz corresponds to the realistic assumption that at the edge of the central calorimeter we will 
have about 2 particles with pt ~ 0.5 to 1 GeV. 

"'^^In this way we essentially fix the value of xp, so that a possible interaction between the upper proton and 
the Pomeron placed in the gap interval (which causes an additional enhanced correction) will not affect the xja 
dependence of R. 




(18) 
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Figure 11: The predictions of the ratio R of (17) for W and T production with (continuous 
curves) and without (dashed curves) enhanced soft rescattering on intermediate partons. 
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Figure 12: The predictions of the ratio R of (17) for the production of a pair of high Et jets with 
(continuous curves) and without (dashed curves) enhanced soft rescattering on intermediate 
partons. 
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that is yS'g^ ~ 1. Moreover, this region corresponds to x ~ 10~^ where the quark densities are 
well known. Therefore inclusive production acts as a good luminosity monitor and the ratio 
R will yield information about the eikonal factor S^. Note that we now probe the survival 
factor at much smaller bt than that corresponding to Figs. 4 and 5. Here we expect the eikonal 
(non-enhanced) survival factor to be S"^ = 0.08. 

On the other hand, diffractive production of a relatively light T is associated with larger 
rapidity intervals available for secondaries, and hence the possibility of more soft rescattering 
with intermediate partons, leads to more enhanced absorption. Indeed, the expected enhanced 
survival factor S"^^ ~ 0.2 — 0.3. For T production the variation of R with ua is weak. The 
smaller rescattering of intermediate partons with y > ua on the lower proton is compensated 
by stronger rescattering of the partons with y < ija on the upper proton. 

Perhaps the most informative probe of 5*^^ is to observe the ratio R for dijet production in 
the region Et 15 — 30 GeV. For example for Et ~ 15 GeV we predict S^^ ~ 0.25, 0.4 and 
0.8 at = —2, and 2 respectively. 

7 Conclusions and Outlook 

Most of the diffractive measurements described above can be performed, without detecting the 
very forward protons, by taking advantage of the relatively low luminosity in the early LHC 
data runs. This allows the use of a veto trigger to select events with no hadronic activity in 
the region corresponding to the large rapidity gap(s). In this way we are able to study central 
exclusive diffractive processes, which should be experimentally accessible at the LHC, that 
probe the various individual components of the formalism used to predict their cross sections. 
The components are sketched in Fig. 1. To summarize, the gap survival factor, S'^, caused by 
eikonal rescattering may be studied as indicated in Figs. 4, 5 and 11, and the possible enhanced, 
contributions as shown in Figs. 9, 11 and 12. The relevant unintegrated gluon distribution, 
fg, can be constrained by observing T production, see Fig. 6, and the QCD radiative effect, T, 
may be checked by observing exclusive two- and three-jet events, see Figs. 8. 

In the first LHC runs it may be difficult to measure the absolute values of the cross sections 
with sufficient accuracy. For instance, it will take time to determine the luminosity with a 
precision better than, say, 10%. Thus the measurements of the ratios of the rate of events with 
and without rapidity gaps (such as (VT-t-gaps/VF inclusive) and {Z+ga,ps/ Z inclusive) etc.) will 
be more reliable in the early data runs. 

When the forward proton detectors are operating, even at moderate integrated luminosity^"^, 
much more can be done. First, it is possible to measure directly the cross section dPaso/dtdMx 
for single diffractive dissociation, pp — > and also the cross section (PcrY^pE/dyidy2 for soft 

central diffractive production, pp ^ p + X + p. These measurements will strongly constrain 

^^That is, before these detectors accumulate enough luminosity to probe for new physics signals. 
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the models used to describe diffractive processes and the effects of soft rescattering. The 
predictions of a very recent detailed model can be found in Figs. 20 and 23 of [47]. It tiirns 
out that (P(JBPE/dyidy2 is particularly sensitive to the detailed partonic content and sizes of 
the various diffractive eigenstates, see Fig. 20(d). 

Next, a study of the transverse momentum distributions of both of the tagged protons, 
and the correlations between their momenta, pa and pt2, is able to scan the proton optical 
density (opacity) [55, 30] (see also [56, 57]). In principle, there are additional smearing effects 
caused by the intrinsic transverse momentum spread of the proton beams, see for instance, [58]. 
Incorporation of these effects requires detailed studies, including, in particular, the detector 
resolution. Therefore the predictions based on measurements of the proton transverse momenta 
should allow for this smearing. 

We emphasize that the selection of central exclusive dijet production in the kinematical 
region corresponding to the sought-after Higgs signal {Ex ~ Mh/2) provides an ideal "standard 
candle" . At leading log accuracy, this process includes all the components of the theoretical 
formalism used to predict the central exclusive Higgs signal; the same parton densities in 
the same kinematical region, the same gap survival factors S'^ (and S^^) and the same QCD 
radiative effect T. The only small difference arises at next-to-leading log order. It is caused by 
the interference between the large-angle soft gluon radiation from the initial active gluons and 
the coloured composition of the dijet system. This effect should be evaluated allowing for the 
detector acceptance and the selection cuts used in the experiment. 

Finally note, that when the statistics allow, a valuable check of the formalism of central 
exclusive diffractive production will come from the central diphoton production, see for instance, 
[59]. The first CDF results [60] are quite encouraging. 
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